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3.3.6,  SORTING  ROUTINE. 

3. 3*6.1.  INTRODUCTION.  The  sorting  routine  sorts  a  list,  stored  as 
array  A,  into  ascending  order. 

3«3,6,2,  INPUTS.  The  sorting  routine  has  to  be  supplied  with  the  appro¬ 
priate  numerical  values  of  the  following  parameters! 

-  beginaddress  of  the  list  Af  i.e.  ii 

-  sndaddress  of  the  list  Af  i.e. 

-  desired  ordering,  i.e.  asoending  or  descending 

-  length  of  the  elements  of  A 

-  the  relative  position  of  the  reference  elements 


in  case  q  has  to  be  adapted! 

-  previous  sorting  time 

-  sign  of  change  of  q  on  previous  run. 

3»3*6.3.  PROCESSING.  The  procedure  sorts  the  elements  of  the  list  A 
[ii  »  JJ]  into  ascending  order.  It  continually  splits  the  list  into 
parts,  such  that  all  elements  of  one  part  are  less  than  all  elements  of 
the  other,  with  a  third  part  in  the  middle  consisting  of  a  single  element. 
An  element  with  value  t  is  chosen  as  pivot  element,  i  and  j  give  the 
lover  and  upper  limits  of  the  segment  being  split.  After  the  split  has 
taken  place  a  value  p  will  have  been  found  such  that  A  [pj  »  t  and 
A  £l  1  4:  t  A  [J  J  for  all  I,  J  such  that  i£  I  ^  p^- J  dr  j. 

The  program  then  performs  operations  on  the  two  segments  A  (  i  ;  p  -  “jj 
and  A[p  +  1  *  jj  as  follows!  the  smaller  segment  is  split  and  the  posi¬ 
tion  of  the  larger  segment  is  stored  in  the  IL  and  IU  arrays  (IL  and  IU 
are  mnemonics  for  Index  Lover  and  Index  Upper).  Their  demension  is  slog 
( j j  -  ii  +  2).  If  the  segment  to  be  split  has  q  or  fever  elements  it  is 
sorted  directly,  instead  and  another  seg~ment  is  obtained  from  the  IL 
and  IU  arrays,  rfhen  no  more  segments  remain,  the  list  A/fiitJJ./  As  com¬ 
pletely  sorted.  The  optimal  value  $  of  the  parameter  q  depends  on  the 
initial  ordering  present  in  the  list  AfiiijjJ  ,  $  being  higher  if  the 
initial  ordering  is  better.  For  random  lists,  a  value  of  6  or  10  for  q 
will  prove  to  be  optimal.  For  non  random  lists,  q  can  be  adapted  by  in¬ 
creasing  it  depending  upon  part  sorting  times. 

The  optimal  value  of  t,  the  value  of  the  pivot  element,  is  the  median  of 
the  list  being  split.  Unfortunately,  this  value  is  known  only  after 
completion  of  the  sorting  process.  Therefore,  an  estimate  of  the  median 
has  to  be  made.  Satisfactory  performance  is  obtained  by  estimating  t 
from  the  median  of  a  sample  of  three  elements  of  A.  These  elements  are 
the  first,  middle  and  last  element  of  the  list  A, 

For  random  lists,  the  method  of  Het  Mathematisoh  Centrum,  which  uses  two 
pivot  elements  instead  of  one,  is  faster,  but  for  non  radom  lists  it  is 
not  faster,  Therefore,  it  is  not  discussed  here.  It  is  a  simple  matter 
to  change  the  sorting  routine  into  one  that  operates  on  equidistant  memory 
elements  instead  of  on  successive  memory  elements.  Also,  the  sorting 
routine  can  easily  be  generalised  into  one  that  sorts  either  into  ascen¬ 
ding  or  descending  order. 
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3. 3*6,3.  (continued) 

In  case  the  sorting  routine  has  to  sort  lists  (tracktables  for  instance) 
according  to  the  numerical  value  of  one  of  their  elements ,  (called  the 
reference  element)  instead  of  sorting  single  elements ,  it  becomes  in¬ 
efficient  to  exchange  two  lists  completely,  when  during  the  sorting  pro¬ 
cess  two  elements  have  to  be  exchanged.  In  this  case,  the  original  lists 
to  be  sorted  are  replaced  by  lists  with  only  two  elements s  one  element 
representing  the  numerical  value  according  to  which  the  sorting  is  to 
take  place,  and  one  representing  the  location  (starting  address)  of  the 
original  lists.  When  the  sorting  of  the  two-element  lists  is  complete, 
the  original  lists  can  be  placed  in  their  proper  place,  if  desired,  thus 
minimising  the  number  of  exchanges. 

The  splitting  of  the  list  into  parts  such  that  all  elements  of  one  part 
are  less  than  all  elements  of  the  other,  proceeds  as  follows:  As  pivot 
element  the  median  value  of  the  first,  middle  and  last  element  is  chosen* 
Further  more,  these  three  elements  are  placed  in  ascending  order. 

Then,  starting  from  the  next-to-last  element  the  list  is  searched  upward 
for  an  element  less  than  or  equal  to  the  pivot  element.  When  such  an 
element  A  J  has  been  found,  the  list  is  searched  downward  starting  froi 
the  second  element  for  an  element  greater  than  or  equal  to  the  pivot 
element.  When  such  an  element  A  [ IcJ  has  been  found,  A  CtJ  and  A C\7  are 
exchanged,  provided  that  K£L,  and  the  search  for  elements  to  be  exchan¬ 
ged  is  repeated  in  the  list  AfK+1  j  L-lJ/  ,  Because  the  middle  element 
is  equal  to  the  pivot  element,  automatic  stopping  is  assured, 

3. 3.6,4.  OUTPUTS. 

-  the  sorting  time 

-  value  of  q 

-  sign  of  change  of  q 

-  beginaddress  of  the  list  A 

-  endaddress  of  the  list  A 

3« 3 .6.5,  SPECIAL  REQUIREMENTS .  When  it  is  required  to  sort  a  list  with 
a  good  initial  ordering,  it  is  advantageous  to  have  an  adaptive  sorting 
routine.  Therefore,  the  sorting  routine  shall  be  adaptive,  i,e.  q  shall 
remain  constant,  or  be  increased  or  decreased  by  one,  depending  upon 
the  previous  and  present  sorting  time,  and  the  sign  of  change  of  q  on 
the  previous  run. 

The  inixal  value  of  q  should  be  chosen  between  6  and  10, 

For  short  lists  (  —  40  elements)  it  iB  felt  that  adaptation  will  not  be 
efficient,  because  the  possible  gain  in  sorting  time  will  be  balanced  by 
the  time  needed  for  adaptation. 
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iUMBadJAia  3ort(a,  11 » j j»q#w)  i 

array  A; 

frfgln  isjlL  *tz9 **»**» y 1 1 9 1 1 }  inf  re  r  Utim  irriY  imu[0tv]i 

switch  8»J**  Ll ,  L2f  L3  y  LA,  L5  | 
ns-  0;  it.  iij  j <-  jjj  goto  LA) 

Lli  pi-  (l4j)  djv  2 )  u«  A[p]j  k: -i  )  L:»  j) 
xi-  A[l ]; 
if  x>t  then 
begin 

yi-t|  A[p]i-ti-x|  A[i]i-xi«y 

end  { 

**“  A[j]| 

IX  z<t  then 
begin 

y*-t>  A[p]«-ti-«i  A,[j]:-zi«y| 

XL  x>t  then 
begin 

yi -t |A[p] i«t i-xj  A[i]s-xi-y 

end 

find.; 

L2:  L:«L-1; 

Zfci  «a[l]| 

If  zz>t  then  goto  L2; 

L3 1  ks -k+l ; 

xxi *A[k  ]| 

If  xx<t  then  goto  L3{ 

If  l^<  L  then 
begin 

yt-zz;  a[lJi-zzi-xx|  A[k]i-xxi-y| 
goto  L2 

en<l ;  If  L-l  >  j-k  then 

begin  XL[m]i-  i{  iu[m J: «  L;  li-k  end 

else 

hoi-i IL.[m ] i-  k;  IL'[«]l-  J|  Ji-  L  |  end  ; 

■  I-  flh»l  ) 

LA  I  J— 1  >  q  then  goto  Ll) 

for  Ll-  1+1  step  1  until  j  to 

begin 

xi^a[l]j 
ki-  L-l ) 
xxi-  a[h]» 
if  xx>x  then 
begin 

L5*  A[k^l]i-xx) 

ki«k-l ) 

XL  *2  l  1/ASfi 

bcx;ln 

xxi  -A[k  J) 

II  xx>x  then  go’.o  L5 


find; 

A[kVl  ]«• 


end 

MU 

ml-  »-l )  ^  ^  Q  then 

begin  Ii«IL[«]T  Jt*XUla]|  goto  L4  end 
end  SORT) 
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